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We have experimentally investigated the fluxon dynamics in a narrow Bi2Sr2CaCu2 08+d stack 
with junction length L ~1.8 fim. As an evidence of high-frequency excitation by a collective cavity 
mode, under an (in-plane) external magnetic field, the current-voltage characteristics show promi- 
nent Fiske steps with the corresponding resonance frequencies of 75-305 GHz. Further study of 
fiux-flow resistance oscillation with various c-axis currents clarifies the correlation with Fiske steps 
by distinguishing two different regions i.e., static flux-fiow region at low bias current level and 
dynamic Fiske step region at high bias current level. 

PACS numbers: 72.30.-|-q, 74.25.Qt, 85.25.Cp 
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In stacked Bi2Sr2CaCu2 08-i-d (BSCCO) intrinsic 
Josephson junctions (IJJs), the mutual interaction be- 
tween junctions is expected since the superconducting 
layers, with the layer thickness d — 0.3 nm, are much 
thinner than the London penetration depth Xl= 170 
nmi. Under the appUed magnetic field parallel to the 
layers, Josephson vortices in a stack form lattice config- 
urations that depend on a constant phase shift between 
neighboring layers, ranging from for the rectangular 
lattice to n for the triangular lattice (see Fig. nja))^*^. 

Under the appropriate conditions, the fluxon lat- 
tice will excite the two-dimensional cavity modes in N 
stacked junctions, leading to the emission of electromag- 
netic wave with characteristic frequency fnm ~ mc„/2L 
(c„ is the phase velocity of electromagnetic wave, L is 
the junction length, m 2, 3,. . . ) denotes i-direction 
mode, and n {—1, 2, 3,. . . N) denotes the stacking direc- 
tion mode)^. The c„ is given by 
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with the Josephson penetration depth Aj, the coupling 
parameter S, the junction number A'' and the Josephson 
plasma frequency WpjiS*^. Among the N different modes 
along the stacking direction, the mode with the lowest 
velocity c^v is usually stimulated by a triangular vortex 
lattice; while for the mode with the highest velocity ci, 
a rectangular lattice is favorableSii. 

Being driven by a c-axis bias current, the moving 
fluxon lattice generates flux-flow resistance (PER) in the 
junctions. Recent experimenta&SiiS and numericaUiii^ 
studies show that the motion of traveling fluxon lattice in 
the layers is reflected into the periodic oscillation of PER 
under the low bias current by considering the dynami- 
cal matching between the vortex lattice and the sample 



edges. The iJo/2-period oscillation of FFR is interpreted 
as a result of the formation of triangular lattice in a long- 
junction stack, where Hq is the field for adding one fiux 
quantum per one junction {Hg = <^q/Ls, $g being the 
flux quantum and s the layer periodicity along the c— axis 
i.e. 1.5 nm)S4ii. Recent studies indicate that when the 
junction size is reduced down to a few ^uv and approaches 
the short-junction stack limit {L < the oscillation 
period becomes predominantly Hq rather than Hq/2 due 
to the deformation of the Josephson vortex lattice by 
strong interaction with junction edge o^i^*^ . The result 
may suggest a possible existence of collective vortex mo- 
tion such as a rectangular vortex lattice in a narrow stack. 
However the recent PPR research is limited at very low 
bias current levels such as less than 1% of critical current 
at zero fleld Icq, and it is still obscure whether the vortex 
dynamics studied by PPR is related with high-frequency 
excitation or not. 

In this paper, stimulated by the above-mentioned re- 
searches, we conflne ourselves to study of the singularities 
in the I — V curves and the correlation with the current- 
dependent PPR oscillation in the narrow BSCCO stacks 
with L ~1.8 /im and large junction number N » 1, 
exploring the high-frequency excitation by a possible col- 
lective fluxon motion. 

In our experiments, a BSCCO whiskeaii with a clean 
and flat surface was fabricated to be an in-line type of 
IJJs stack with the length of 1.8 /im by a focused ion 
beam. A schematic diagram of the stack is shown in 
the inset of Fig. ^b). Here, L and D denote the junc- 
tion length perpendicular to the magnetic fleld and the 
depth parallel to the field. The fabricated IJJs stack had 
a thickness of about 180 nm, i.e., containing about 120 
junctions. The values of critical temperature (Tc) was 81 




vortex lattice in a stacked BSCCO IJJs. The trianeular lattice 1 2 3 4 5 
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corresponds to the out-of-phase case, while the rectangular 
lattice corresponds to the in-phase situation, (b) Scanning 
ion-beam microscope image of BSCCO stack fabricated by 
the focused ion beam etching. Schematic diagram of the stack 
is shown in the inset. Magnetic field H was applied along the 

longer side (D) of junction to enhance edge effect. 1 2 3 4 5 

Normalized field h 
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K, and the values of Ico were about 350 /uA at 10 K and 
about 195 iiA at 50 K respectively. The electric transport 
properties were measured with a four-terminal configu- 
ration using a Physical Property Measurement System 
(PPMS, Quantum Design), which can supply magnetic 
field up to 9 T. In order to enhance the edge effect of 
the sample on moving Josephson vortices, we applied the 
magnetic field parallel to the a&-plane along the longer 
side (D) of the BSCCO stack (see Fig-db)). The sam- 
ple was mounted on a rotatable holder with a resolution 
better than 0.005°. The m-plane alignment was precisely 
adjusted by the angular dependence of FFR under exter- 
nal magnetic field. 

The Josephson penetration depth Aj is given by Aj ~ 
[$o/27r/^oic(ie// + 2A^/de//)]^^^ with the eflfective val- 
ues deff = Asinh(c?/A), t^fj = t + 2A tanh(d/2A), 
the (in-plane) magnetic penetration depth A(r) = 
A,6(0)/v/l - (r/Tc)4, by assuming Aa6(0)=170 nm, the 
thickness of superconducting layers and insulating lay- 
ers c?=0.3 nm and t=1.2 nmi^. The Aj of our sample 
was calculated as 0.27 fim at 10 K (0.34 /zm at 50 K). 
The sample width of 1.8 /Ltm is about 6.6 times (about 5.3 
times at 50 K) as large as the calculated Aj at 10 K, thus 
it can be still regarded as a long-junction stack; however, 
as we applied magnetic fields perpendicular to the nar- 
rower side, which is different from the conventional case, 
we use the term 'narrow' to specify our sample. 

Figure|2ta) displays the I—V characteristics under var- 
ious magnetic fields parallel to the layers at 50 K. In or- 



FIG. 2: (Color online) (a) I ~ V characteristics at various 
fields from 0.765 T (/i = f) to 3.825 T {h = 5) with an 
interval of 0.0765 T (O.f h) at 50 K. Bold lines show I -V 
curves at fields h= integer and half integer where Fiske steps 
are clearly observed, (b) The magnetic field dependence of 
the current amplitudes of first step and second step at 50 K. 

der to indicate a periodic modulation of the I ~V curves 
with field, we normalized the magnetic fields in graphs 
of Fig. 2 by the period Hp— 0.765 T, which corresponds 
to Hq calculated for this sample. Then we can regard 
the normalized field h = H/Hp{= H/Hq) as Josephson 
fiuxon number per unit junction. The magnetic field was 
applied in the range from 0.765 T [h = 1) to 3.825 T 
(h = 5) with an interval of 0.0765 T (O.l/i). 

With increasing magnetic field, we found clear current 
steps that were regularly developed in the I—V curves, as 
shown in Fig.l^^a). These steps were identified as Fiske 
steps that had been observed as a strong enhancement of 
superconducting current when the Josephson frequency 
{uj = 2eV/H) matches the resonant frequency of elec- 
tromagnetic cavity modes excited in junctionsi^iiSiiLi^. 
The asymptotic voltage positions of step series are given 

by 

V = m7V($oc„/2L). (2) 

Using the voltage position of the first Fiske step (18.7 
mV) defined by local maximum dl/dV, the characteristic 
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velocity was estimated to be c„ = 2.71 x 10^ m/sec using 
Eq. lO with m=l, iV=120, and L=1.8 fiiR. Since the 
series of Fiske steps are observable up to the 4th order in 
Fig. |2Ia), the corresponding resonance frequency lies in 
the range of 75-305 GHz. 

According to Eq. using a set of the experimental 
parameters (the current density jc=l-015 kA/cm^ at 50 
K, the thickness of insulating layer t—1.2 nm, the cou- 
pling constant 5* « 0.5, and Aj=0.76 /im), the velocities 
of the lowest mode and the highest mode are calculated 
as ci2o = 2.62 x 10^ m/sec and ci = 1.86 x lO"^ m/sec, 
respectively. Our result, estimated from the Fiske step, 
is comparable with the velocity of the lowest collective 
cavity resonance mode, in agreement with Refs. ^3 and 

E 

The observed odd and even steps obviously have differ- 
ent dependences on the magnetic field. Fig. |2Ib) shows 
the field dependence of the step height (SIps) for the first 
and second steps. With increasing field, each height of 
the first and the second order steps oscillated with the 
same period Hq- However, the maxima of the first order 
step appeared at /i=integer, and the minima at h—hali- 
integer; while the second order step exhibited opposite 
behavior. Note that this is very similar to the behavior 
of a single junctioni^ti^ except that large current-step 
voltage due to the contribution from all junctions. This 
is very important because such a single junction behav- 
ior is possible when all junctions are evenly excited by a 
collective cavity resonance in stacked junctions. 

For further understanding of Fiske steps, we also stud- 
ied the correlation between Fiske steps and the FFR os- 
cillation as function of magnetic field, which has been 
regarded as a powerful tool for this purposa^iSiiSiiLiSi. 

At extremely low c-axis current I — 1 /lA (0.52 % of 
Ico), as shown in Fig.|31 there are two oscillation regimes 
of FFR: one is below h ~ 2.5 and the other is over h — 2.5 
with the oscillation period of Hq (see Part A and B in 
Fig. |3J|. The details for the change in oscillation period 
by magnetic field will be elucidated with the model based 
on the edge current elsewhereiS. 

Having noticed in Fig. |2Ia), there is not any step at 
such a low current level, we intended to measure the FFR 
at much high current levels. Surprisingly, we have ob- 
served a new anomaly of FFR oscillation at higher bias 
currents as shown in Fig. In particular, for high field 
h > 2.5 (Part B in Fig. ^ it is found that the peaks 
and the local minima of FFR oscillation were inverted 
by varying bias current as marked by the arrows, with 
the oscillation period Hq unchanged. Note that the in- 
version from the peaks to the local minima of FFR oscil- 
lation implies a change of the matching between the edge 
and Josephson vortices as well as an enhanced conduc- 
tance. As is well known, such conductance enhancement 
in a Josephson junction, usually resulting from resonance 
in the junction, should manifest itself as a current steps 
in the I — V characteristics. 

To confirm this, we compared the current-dependent 
anomaly of FFR oscillation in Fig.Owith the I—V curves 
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FIG. 3; (Color online) FFR anomalies due to c-axis bias cur- 
rent of 1-10 fiA at 50 K, where each curve is shifted by 1 kQ 
for clarity. The magnetic fields are normalized by Hp = 0.765 
T. With increasing current, there are two kinds of features; 
Part A: two small peaks transform to one peak with period 
Ho at h < 2.5. Part B: the inversion from peaks to local 
minimums at ft > 2.5. The dash dot line in Part B shows the 
boundary of two different regimes with c-axis bias current. 

in Fig.[21[a). Looking at the FFR curves marked by the 
arrows at the normalized field h — A and 5 in Fig. |2| the 
clear inversion from peaks to local minima took place 
when the bias current was varied from 5 //A to 9 fj,A. 
Considering the corresponding curves in Fig. Ufa) with 
the same field and bias current range, we see two Fiske 
steps appearing at the same voltage in the I — V curves, 
as marked by the arrows at the fields h — A and 5. As 
one can see in Fig. [SJ there are several more inversion 
points at different normalized fields, and the correspond- 
ing steps in Fig.[2Ia). 

Different from the high field regime h > 2.5 where the 
FFR oscillates with the period Hq, in the low field regime 
h < 2.5 where the FFR oscillation period is smaller than 
Hq and comparable with Hq/2, there is neither inversion 
of the FFR oscillation nor steps up to 10 /j,A. In the case 
of the field of /i = 2 in Fig. |2ta), the first Fiske step 
appears at the current around 20 /xA. The experimental 
results show the inversion of FFR oscillation near 20 fiA 
after FFR becomes i/o-oscillation, i.e., two peaks merge 
into one peak (see the arrows in Part A of Fig. |2J|. Ac- 
cordingly, it is clear that the Fiske steps appear only in 
the regime that FFR oscillates with a period of Hq and 
that there is the inversion of FFR oscillation from peaks 
to local minima. 

Noticeably, after the inversion of FFR oscillation, there 
is distortion of i?o-oscillation with /=8 /xA around h—A 
and ft,=4.5 where the first and the second steps appear in 
I — V curves. This can be understood with the following 
explanation. When the first and the second Fiske steps 
coexist at same current level, FFR will reflect the behav- 
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ior of these Fiske steps under magnetic field. Then the 
product of two (odd and even step's) field dependence 
factors (refer to Fig. [21 (b)) should consequently result 
in i7o/2-period oscillation of FFR. This is in fact the 
case at which Ustinov and Pcdcrscn could observe Hq/2- 
oscillation in their simulation for a single long-junction^S. 
However in our experiments, some steps don't coexist at 
the same current due to large voltage intervals between 
Fiske steps and the reduced step amplitude at 50 K. 
Therefore the FFR oscillation was mainly affected from 
one step (odd or even), with the dominant Tfo-period at 
the current near Fiske steps. Nevertheless, the neighbor- 
ing steps caused observable distortion of FFR. 

Having understood the correlation between the FFR 
oscillation and the I — V characteristics, we can distin- 
guish between the two different regimes divided by dash 
dot line in Part B of Fig. O In the low bias region be- 
fore the occurrence of inversion of FFR, i.e., before the 
appearance of Fiske steps, the Josephson fluxon dynam- 
ics is mainly determined by the edge pinning effect and 
the coupling interaction between the layersiiii^. There- 
fore the information of ffuxon lattice can be probed by 
the FFR oscillation measurement to some extent. As a 
matter of fact, such oscillation of FFR in this regime (es- 
pecially at the low bias current and high field) is just 
same as Ic modulation with magnetic field. 

On the other hand, in the inversion region of FFR os- 
cillation accompanying the Fiske steps at relatively high 
bias current range, the dynamic resonant fiuxon motion 
is dominant due to the interaction between the travel- 
ing fluxon and the cavity mode excitation. Thus in the 
high bias region showing the Fiske steps, the fiuxon lat- 



tice can't be determined simply by the measurement of 
FFR oscillation. We can only try to figure out the lat- 
tice structure by comparing the characteristic velocity 
obtained from the Fiske steps measurement and Eq. 

With low-temperature scanning electron microscopy 
(LTSEM), Quenter et ali^ also observed similar different 
regimes . For the bias points close to origin (in a finite 
voltage state without Fiske steps), the LTSEM image 
showed the static distribution of the Josephson current 
in the presence of an external magnetic field. For the high 
bias point close to a voltage of Fiske resonance, a clear 
standing-wave pattern is observed due to the superposi- 
tion of the traveling wave and the refiected waves. These 
are consistent with our results. 

In summary, the current-voltage characteristics of 
BSCCO IJJs stacks with L ~1.8 fim under an exter- 
nal magnetic field showed pronounced Fiske steps as an 
evidence of high-frequency excitation. The alternative 
appearance of even and odd Fiske steps resembles the be- 
havior of a single junction, although there are more than 
100 junctions in one stack. The observed Fiske steps and 
their mode of the collective cavity resonance suggest that 
all junctions in such a narrow stack can be synchronized 
by the Fiske resonance. Further measurements on the 
field-dependent FFR with various c-axis currents clar- 
ify the correlation with Fiske steps by distinguishing two 
different regions i.e., static flux-flow region at low bias 
current level and dynamic Fiske step region at high bias 
current level. 
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